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Periodically driven quantum systems provide a novel and versatile platform for realizing topolog¬ 
ical phenomena. Among these are analogs of topological insulators and superconductors, attainable 
in static systems; however, some of these phenomena are unique to the periodically driven case. 
Here, we show that disordered, periodically driven systems admit an “anomalous” two dimensional 
phase, whose quasi-energy spectrum consists of chiral edge modes that coexist with a fully lo¬ 
calized bulk - an impossibility for static Hamiltonians. This unique situation serves as the basis 
for a new topologically-protected non-equilibrium transport phenomenon: quantized non-adiabatic 
charge pumping. We identify the bulk topological invariant that characterizes the new phase (which 
we call the “anomalous Floquet Anderson Insulator”, or AFAI). We provide explicit models which 
constitute a proof of principle for the existence of the new phase. Finally, we present evidence 
that the disorder-driven transition from the AFAI to a trivial, fully localized phase is in the same 
universality class as the quantum Hall plateau transition. 


I. INTRODUCTION 


Time-dependent driving opens many new routes for 
realizing and studying topological phenomena in many- 
body quantum systems. Recently, an intense wave of 
activity has developed around exploring the possibilities 
of using periodic driving to realize “Floquet topological 
insulators”, i.e., driven system analogues of topological 
insulators [1-19] in a variety of solid state [20], atomic, 
and optical contexts [21, 22]. Beyond these analogies, 
driven systems may also host their own unique types of 
robust topological phenomena, which have no analogues 
in non-driven systems [3, 23-27]. The latter will be at 
the heart of the present work. 

In static (time-independent) two dimensional systems, 
the appearance of chiral edge states is intimately tied to 
the topological structure of the bulk Bloch bands, cap¬ 
tured by the so-called Chern number [28]. This well- 
known bulk-edge correspondence breaks down in period¬ 
ically driven crystalline systems, where chiral edge states 
can exist even if the Chern numbers of all the bulk bands 
are zero [3, 24]. Such anomalous edge states are captured 
instead by a topological invariant that characterizes the 
time evolution operator of the bulk wave functions [24]. 
A system exhibiting this anomalous behavior has been 
recently realized using microwave photonic networks [29]. 

This unique topological phenomenon opens new possi¬ 
bilities, which are inaccessible in static systems. For ex¬ 
ample, it is well-known that Bloch bands with non-zero 
Chern numbers cannot be spanned by a complete basis 
of localized Wannier functions [30, 31]. Correspondingly, 
when quenched disorder is introduced, not all the states 
in a Chern band can be localized; delocalized states must 
exist at least at one value of energy in the band [32]. In¬ 
tuitively, this can be understood from the fact that the 
chiral edge states in the bulk band cannot “terminate” 
without hybridizing with a delocalized bulk state. In 


FIG. 1. The anomalous Floquet-Anderson insulator (AFAI), 
in a disordered two-dimensional periodically-driven system 
with time-dependent Hamiltonian H(t). In the AFAI phase 
all bulk states are localized, yet the system hosts chiral prop¬ 
agating edge states at all quasienergies. The nontrivial topol¬ 
ogy of the phase is characterized by a nonzero value of the 
winding number defined in Eq. (2) 


contrast, a periodically driven system can exhibit chiral 
edge states even when all the Chern numbers are zero. 
Moreover, due to the periodicity of quasienergy, it is in 
principle possible for a chiral edge state to wrap around 
the entire quasi-energy zone without terminating at a de¬ 
localized bulk state. This leads us to hypothesize that, 
in a disordered periodically-driven system, robust chiral 
edge states may coexist with an entirely localized bulk. 
In such a system, which we term an anomalous Floquet- 
Anderson insulator (AFAI), the chiral edge states form 
a uni-directional one dimensional system, whose dynam¬ 
ics is decoupled from the bulk at all quasienergies. This 
situation defies the standard intution from strictly one¬ 
dimensional systems that must have an equal number of 
right and left moving modes. But can an AFAI state re¬ 
ally exist? And if so, what are its physical consequences? 

In this work, we explore the AFAI phase in period¬ 
ically driven, two dimensional disordered systems. We 
construct explicit models that demonstrate its existence, 
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and discuss its topological characterization and its phys¬ 
ical properties. Strikingly, the AFAI hosts a unique non- 
equilbrium topological transport phenomenon: quantized 
charge pumping in a non-adiabatic setting. Essentially, 
if all the states in the vicinity of the edge are occupied 
by fermions (to a distance of several times the bulk lo¬ 
calization length), the uni-directional edge states carry a 
current whose long-time average is quantized in units of 
one particle per driving period. Importantly, disorder is 
essential for the quantization of pumping in the AFAI; 
absent the disorder, generically there is no quantization 
due to the presence of delocalized states in the bulk. 

Quantized pumping is well-known from the work of 
Thouless on adiabatically driven one-dimensional sys¬ 
tems. however, unlike in the Thouless pump, in the 
AFAI, the driving frequency is not required to be small in 
order to observe the quantization of the current. The rea¬ 
son the current at the edge of an AFAI can remain quan¬ 
tized even when the adiabatic condition is violated is that 
the two counter-propagating edge modes that carry the 
current are spatially separated [33]; hence, they cannot 
backscatter into each other even if the driving frequency 
is not small. 

This paper is organized as follows. In Sec. II we intro¬ 
duce the defining properties of the AFAI; we discuss the 
topological invariant characterizing the AFAI, and show 
that the AFAI exhibits edge modes at every quasi-energy. 
In Sec. Ill we show how the edge mode structure leads 
to quantized charge pumping. We then demonstrate, in 
Sec. IV, the appearance and robustness of an AFAI in a 
simple, tractable model. In Sec. V we conduct a numeri¬ 
cal study of a wider class of models exhibiting the AFAI 
phase. We numerically demonstrate the properties dis¬ 
cussed in sections II and III. At strong disorder, we find 
a topological transition between the AFAI and a “trivial” 
Floquet insulator where all states are localized (includ¬ 
ing at the edges); we speculate that the transition is in 
the same universality class as the quantum Hall plateau 
transition, and corroborate this using our numerical re¬ 
sults. 


II. THE AFAI: TOPOLOGICAL INVARIANTS 
AND EDGE STATES 

We begin by defining the AFAI, and introducing the 
topological invariant which characterizes it. The defin¬ 
ing characteristic of the AFAI phase is the peculiar re¬ 
lationship between its bulk and edge mode spectra: in 
the AFAI phase all bulk Floquet states are localized, 
yet the system still hosts topologically-protected chiral 
modes along its edges. 

We consider a two-dimensional system of non¬ 
interacting particles with a time-periodic Hamiltonian, 
H{t) = H(t + T), where T is the driving period. No spa¬ 
tial translational symmetry is assumed. The interesting 
aspects of the AFAI phase are revealed by comparing the 
Floquet operators U(T) = Te~ l ^ o dtH( d ) for toroidal and 


cylindrical geometries. In an AFAI phase, all the eigen¬ 
states of U (T) on a torus are localized. However, as we 
will show, in a cylindrical geometry there are eigenstates 
which are localized at the boundaries of the cylinder, but 
delocalized along the direction of the boundary, at every 
quasi energy 0 < e < ft = 2 tt/T. 

The topological invariant which describes the AFAI is 
a generalization of the “winding number” introduced in 
Ref. [24]. As a first step in constructing the topological 
invariant, we define an associated, time-periodic evolu¬ 
tion operator for the system on a torus: 

U £ (t) = U(t) exp (iH^t) , (1) 

with ^ log U(T). Note that, by construction, 

U £ (T) = 1. The explicit dependence on 5 in the above 
definitions comes from the necessary choice of a branch 
cut for log; we use a definition such that —i \oge lx = x 
if X G [0 ,eT) and —iloge* x = x — 27r if x F [eT, 2 tt). 
As an additional ingredient, we also consider a family of 
time-dependent Hamiltonians #(©,£) and the associated 
evolution operators Z^(@,t), in which constant (time in¬ 
dependent) fluxes 0 = ( 0 xi 0y ) are threaded through the 
torus [34]. 

With these definitions at hand, we can define the 
“winding number” 

r T /* d 2 0 

W e = J dt J 8^ ^ dtUe H^UeMldeM) • 

( 2 ) 

The winding number W £ is an integer, which can in prin¬ 
ciple depend on the quasi-energy 6. Note that in order for 
W £ to be well defined, the quasi-energy 5 has to remain in 
a spectral gap of 1/(0, T) for every value of the threaded 
fluxes 0 (otherwise, the operator U £ is discontinuous as 
a function of 0). We argue that for a large enough sys¬ 
tem, almost all values of 5 satisfy this requirement. This 
is because, upon changing fluxes 0 X and 6 y , the quasi¬ 
energies of the localized bulk states only change by an 
amount proportional to where £ is the localiza¬ 

tion length and L is the linear system size. In contrast, 
the average level spacing is proportional to 1/L 2 . 

Next, we show that if all the eigenstates of U(T) are 
localized, then the invariant W £ is in fact independent 
of 6. This follows from the relation between the winding 
number W £ and the Chern numbers characterizing the 
eigenstates of U(T) [24], 

W £1 -W £2 =C £U£2 . (3) 


In the above, C £l ^ £2 is the total Chern number of the 
eigenstates with quasi-energies between e 1 and £ 2 : 
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where Pq 1,£2 ^ is a projector onto the eigenstates of 
U (0, T) with quasi-energies between £1 and £ 2 - If all the 
bulk eigenstates are localized, C £l:£2 = 0 [35]. There¬ 
fore in this case by Eq. (3), W e% = W £2 for every pair of 
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FIG. 2. Edge states and spectral flow in the AFAI. a) The 
parent phase of the AFAI is a clean system without disorder, 
where all Floquet bands have Chern number zero but the 
winding number (2) is non-zero in all gaps. In a cylinder 
geometry, chiral edge state propagate along the upper and 
lower boundaries, only at quasi-energies within the bulk gaps, 
b) The corresponding spectrum, shown as a function of the 
conserved circumferential crystal momentum component, c) 
When disorder is added, all bulk states become localized while 
the chiral edge modes on the cylinder persist. When all states 
are filled near one end of the cylinder, a quantized current 
flows along the edge, d) With disorder, crystal momentum is 
no longer a good quantum number. However, the spectrum of 
states localized near the upper edge, displayed as a function 
of the flux 0 X threaded through the cylinder, clearly displays a 
non-trivial spectral flow. The spectral flow fully winds around 
the quasi-energy zone, accounting for the quantized pumping 
in the AFAI phase. 


quasi-energies. Then, we can drop the subscript e, and 
refer to the winding number simply as W. 

We thus define the AFAI as a time-periodic, disordered 
system in which (1) all the bulk Floquet eigenstates are 
localized, (2) the quasi-energy independent winding num¬ 
ber W is non-zero. Below, we argue that the boundaries 
of the AFAI necessarily support chiral edge states at ev¬ 
ery value of the quasi-energy. 

Having defined the AFAI phase, we briefly mention 
how it can be reached. A good starting point for ac¬ 
cessing the AFAI is a clean (translationally invariant) 
Floquet-Bloch system, for which all the Chern numbers 
of U(T) vanish, but with W £ ^ 0 for any e within one 
of the gaps in the quasi-energy spectrum. A schematic 
quasi-energy spectrum of such a system in a cylindrical 
geometry is shown in Fig. 2(b). We then add a static, 
spatially disordered potential to H(t). We will argue be¬ 
low that all the bulk states are generically localized even 
for arbitrarily weak disorder. The crucial point is that 


the winding number W need not be zero, even if all the 
bulk states are localized. A specific solvable model, which 
serves as a proof of principle for the existence of the AFAI 
phase, is given in Sec. IV. 

We finish this section with a discussion of the edge 
structure of the AFAI. In the clean limit, there are chiral 
edge states in any bulk quasi-energy gap with a non-zero 
winding number [24]. Clearly, these edge states cannot 
localize when disorder is added. Moreover, intuitively, 
if all the bulk states are localized, the chiral edge states 
must persist even within the bulk bands. To see this, con¬ 
sider a system in a cylindrical geometry. Upon inserting 
a flux quantum through the hole of the cylinder, the chi¬ 
ral edge states exhibit a non-trivial “spectral flow”: i.e., 
even though the spectrum as a whole is periodic as a 
function of flux, every state evolves into the next state in 
the spectrum [Fig. 2(d)]. The spectral flow cannot ter¬ 
minate in the bulk bands. Since all the bulk states are 
localized, they are insensitive to the flux, and hence there 
must exist a delocalized, chiral edge state at every quasi¬ 
energy within the bulk bands that “carry” the spectral 
flow. 

To make this argument more precise, we define a topo¬ 
logical invariant that directly characterizes the spectral 
flow of the edge states. This topological invariant turns 
out to be equal to the bulk invariant W; we will show 
this in detail in Sec. Ill and Appendix C, where we 
demonstrate that both invariants are related to quan¬ 
tized charge pumping along the edge. 

To construct the edge topological invariant, we con¬ 
sider a cylinder that extends from y = 0 to y = L y , with 
a flux 6 X inserted through the hole of the cylinder. The 
evolution operator of the system on the cylinder is de¬ 
noted by U(t). We now isolate the topological features 
of the edge states by deforming the evolution operator 
in the regions away from the edges, such that the evo¬ 
lution in the bulk takes a simple universal form, while 
the evolution near the edges is unaffected. In particu¬ 
lar, we flatten the bulk evolution such that all Floquet 
eigenstates localized sufficiently far (at least a distance 
^o) from the edges have quasienergy 6 = 0. The resulting 
evolution operator interpolates smoothly between U(t) 
in the vicinity of the edge and U £ (t) of Eq. (1) in the 
bulk (an explicit formulation of the deformation proce¬ 
dure appears in Appendix A). The deformed evolution 
operator takes a block-diagonal form: 

(Ui(T) 0 0 \ 

Ue(T) =010, (5) 

V 0 0 U 2 (T)J 

where in the above, the sub-blocks IA\{T) and U 2 (T) cor¬ 
respond to sites with 0 < y < £q and L y — to < V < L y , 
respectively; the unity block acts on sites with < y < 
Ly — fo. The precise value of fo is not important, as long 
as it is much larger than the bulk localization length of 
the original evolution operator, U (T). The integer-valued 
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“edge winding number” is defined as 


is given by 


^edge 



dO x 

27r 


Tr 


U x {Tj*deMT) 



(6) 


where the sum in the second line runs over all the eigen¬ 
states of Wi(T), and ej are their corresponding eigen¬ 
values. The edge winding number (6) counts how many 
times the spectrum oiU\(T) “wraps” around the quasi¬ 
energy zone, 5 G [0,27r/T), as 0 X varies from 0 to 2tt. A 
schematic example of a spectrum with a non-zero winding 
number is shown in Fig. 2(d). Note that the total wind¬ 


ing number of the system, 


r 2n do* 

JO 2tt 


Tr 


U^deMT) 


must vanish [3]. Hence, the winding numbers of U\{T) 
and U 2 (T) must sum to zero. 

A non-zero n e d ge necessarily implies that there are de¬ 
localized states along the edge; if all states were local¬ 
ized, their quasi-energies would be almost insensitive to 
Q x , and hence n e dge would be zero. Note also that, since 
in the AFAI all the bulk states are localized, changing 
would not change n e d ge ; this amounts to adding a few 
localized states to the spectrum of Wi(T), and cannot 
change its winding number. 


III. QUANTIZED CHARGE PUMPING 


We now discuss the physical implications of the AFAI 
phase. Consider an AFAI placed in a cylindrical geome¬ 
try, as in Fig 2(c). Fermions are loaded into the system 
such that in the initial state all the lattice sites are filled 
up to a distance of £ £ from one edge of the cylinder, 

and all the other sites are empty. Below we show that 
in the thermodynamic limit, the current across a verti¬ 
cal cut through the cylinder, averaged over many driving 
periods, is equal to n e d ge , Eq. (6), divided by the driv¬ 
ing period T. The exact form in which we terminate 
the filled region will not matter, as long as all the sites 
near one edge are filled, and all the sites near the other 
edge are empty. The system thus serves as a quantized 
charge pump, but unlike the quantized pump introduced 
by Thouless [33], there is no requirement for adiabaticity. 

In Appendix C 2 we furthermore show by a direct eval¬ 
uation that the long-time average of the pumped charge 
per driving period is also equal to IF, the bulk invariant. 
In particular, this implies that IF = n e d ge - 

To set up the calculation of the charge pumping in the 
AFAI, we choose coordinates such that x is the direction 
along the edges of the cylinder, and y is the transverse 
direction. We denote the initial many-body (Slater de¬ 
terminant) state, in which all sites up to a distance of 
£ from the edge are filled, by |F(0)). Then, the charge 
pumped across the line x = xq between t — 0 and t = r 


(Q)r = 



dH{9 x ,t) 

d6 x 



(7) 


Here, 0 X is the flux through the cylinder and H(Q x ,t) 
is the corresponding Hamiltonian. For Eq. (7), we use 
a gauge such that on the lattice, every hopping matrix 
element that crosses the line x = xq has a phase of e l9x . 

The initial state |F(0)) clearly does not return to itself 
after a single driving period. Therefore, we cannot expect 
that the pumped current to be identical between different 
periods along the evolution, nor can we expect it to be 
exactly quantized. However, we find that the average 
pumped charge over N periods approaches a quantized 
value Qoq in the limit of a large number of periods, where 
the correction to the quantized value decays as 1/N: 


{Q)nt 
N 


Qoo + O 



( 8 ) 


Here, = W = n e d ge (where IF, n e d ge are the bulk 
and edge topological invariants, respectively, defined in 
Sec. II). Note that Qoo is independent of xo, i.e., the 
charge pumped across any line parallel to the y axis leads 
to the same Qoq. 

In order to compute the charge pumped per period, it is 
useful to express |F(£)) as a superposition of the Floquet 
eigenstates. As we show in Appendix C 1, when averag¬ 
ing the pumped charge over N periods, the contribution 
of the off-diagonal terms between different Floquet eigen¬ 
states decays at least as fast as 1/N. The diagonal terms 
yield a contribution that depends on the evolution over 
a single period, giving 

Qoo = £ rij / dt{ipj (t) | | ipj (t)). (9) 

In the above, \ipj(t)} are the single particle Floquet 
states, which evolve in time as | 'ipjit)} = e~ l£t \4>j(t)) 
(where \4>j(t)) is periodic in time), and nj are the Flo¬ 
quet state occupation numbers in the initial state, nj = 
(T(O)|' 0 ]V ; jl^ r (O ))5 where is the creation operator cor¬ 
responding to |^(0)). [Note that if fermions were ini¬ 
tialized in the Floquet eigenstates, such that nj = 0 or 
nj = 1, we would obtain (Q) nt /N = Q^, without the 
correction terms in Eq. (8)]. 

Straightforward manipulations yield Q^ = 
T JF njdsj / d0 x . At this point, the average cur¬ 
rent per period depends on 0 X . In the thermodynamic 
limit, we expect this dependence to disappear. As in the 
case of the quantization of the Hall conductance [36], we 
average over Q x [37]. We therefore get 

_ r 27 T r\ 

Q °° = ^12 J 0 de ^df x (10) 

Equation (10) relates the average current in a period 
to the spectral flow of the Floquet spectrum as the flux 
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O x is threaded. It is reminiscent of the expression for the 
edge topological invariant, n e d ge , Eq. (6), defined in terms 
of the “deformed” evolution operator U e (T). Below, we 
give a heuristic argument that indeed Q oo = n e d ge , up 
to corrections that are exponentially small in i. A more 
rigorous (but technically cumbersome) derivation of the 
relation between the pumped charge and the bulk invari¬ 
ant is presented in Appendix C. Numerical evidence for 
the quantization of the pumped charge is shown in Sec. V. 

Our strategy in analyzing the pumped charge is to de¬ 
form the evolution operator into the “ideal” form, U e (T) 
of Eq. (5), for which the pumped charge is exactly quan¬ 
tized, and to put bounds on the correction to the pumped 
charge due to the deformation. We define the deforma¬ 
tion process according to Appendix A, with Iq, the width 
of the strip beyond which the quasi-energy spectrum be¬ 
comes flat, chosen such that i ^ £q. Clearly, for the 
deformed evolution operator, rij = 1 for every eigenstate 
ofWi. Therefore, the deformed evolution operator has an 
exactly quantized pumped charge, equal to n e d ge . 

Now, consider the pumped charge of the original (un¬ 
deformed) evolution. We can roughly divide the Floquet 
states that contribute to Eq. (10) into three categories: 

1. States that are localized far from occupied region, 
y^i. For these states, rij is exponentially small, 
and hence their contribution to Qoq is negligible. 

2. States that are localized near the edge, y t. 
These states have rij « 1. Their wavefunctions 
and quasi-energies, and hence their contribution to 
Q oo, are essentially unaffected by the deformation 
process. 

3. States that are localized near the boundary be¬ 
tween occupied and unoccupied sites, y ~ For 
such states, rij is neither close to 0 nor to 1; how¬ 
ever, these states are localized in the x direction 
(as are all the bulk states in the AFAI). Therefore, 
dsj/d0 x of these states is exponentially small, and 
they contribute negligibly to Qoo- 

As 0 X varies, there are avoided crossings in the spec¬ 
trum, in which the character of the eigenstates changes. 
E.g., an eigenstate localized around y\^i may undergo 
an avoided crossing with an eigenstate localized around 
2/2 ^ A When 0 X is tuned to such degeneracy points, the 
two eigenstates hybridize strongly, and do not fall into 
either of the categories discussed above. Such resonances 
affect both d£j/dO x and the occupations rij of the res¬ 
onant states. However, since the eigenstates that cross 
are localized in distant spatial areas, the matrix element 
that couples them is exponentially small. Therefore sig¬ 
nificant hybridization requires their energies to be tuned 
into resonance with exponential accuracy, limiting the 
regions of deviation to exponentially small ranges of 6 XJ 
of order e~ £ ^. The number of such resonances increases 
only polynomially with the size of the system, and there¬ 
fore for L y i £ and L x oc L y , their effect on Qqq is 
exponentially small. 
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FIG. 3. Simple model for achieving the anomalous Floquet- 
Anderson phase. The Hamiltonian is piecewise constant, de¬ 
fined in five equal length segments of duration T/5. During 
steps 1-4, nearest-neighbor hopping is applied along the col¬ 
ored bonds as shown. The hopping strength J is chosen such 
that a particle hops between adjacent sites with probability 
one during each step. In the fifth step, all hopping is turned 
off and a random disorder potential is applied (the same po¬ 
tential is used for all subsequent driving cycles). 


We conclude that, in the thermodynamic limit, all the 
contributions to Qoo in Eq. (10) that are not exponen¬ 
tially suppressed are also exponentially insensitive to the 
deformation process. Therefore, Q^ = n e d ge - 


IV. MODEL FOR AN ANOMALOUS 
FLOQUET-ANDERSON PHASE 

In this section, we study a simple model which allows 
us to explicitly demonstrate the existence and robustness 
of the AFAI phase. We start from a solvable model in¬ 
troduced in Ref. [24], which exhibits perfectly flat bulk 
Floquet bands, and hosts chiral edge modes at its bound¬ 
aries. Adding a specific kind of disorder to this model 
results in localization of all the bulk states, while pre¬ 
serving the edge states; the system is thus in the AFAI 
phase. We then argue that this phase is robust to generic 
small perturbations (i.e., the bulk states remain localized, 
and the chiral edge states persist). 

We consider a system on a square lattice with a 
periodic, piecewise-constant Hamiltonian of the form: 

ffclean© = H n , for < t < *£, n = 

The square lattice is divided into two sublattices, A and 
B (shown as filled and empty circles in Fig. 3, respec¬ 
tively). During each of the first four segments of the driv¬ 
ing, n = 1,..., 4, hopping matrix elements of strength J 
between the A and B sublattices are turned on and off 
in a cyclic, clockwise fashion, as shown in Fig. 3: during 
segment n — 1, 2, 3, or 4, each site in the A sublattice is 
connected by hopping to the site above, to the right, be¬ 
low, or to the left of it, respectively. In the fifth segment 
of the period, all the hoppings are set to zero, and an on¬ 
site potential Sa,b is applied on the A and B sublattice 
sites, respectively. 

We choose the hopping strength J such that ^ = f • 
For this value of J, during each hopping segment of the 
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driving period a particle that starts on one of the sites 
hops to the neighboring site with unit probability. The 
on-site potential, applied only while all hopping matrix 
elements are turned off, is chosen to be Sa,b = 

With this time-dependent Hamiltonian, it is easy to find 
the Floquet eigenstates and quasi-energies. The bulk 
spectrum consists of two flat Floquet bands with quasi¬ 
energies =b^p, with the corresponding eigenstates local¬ 
ized on either the A or B sublattice. The winding number 
invariant can be computed for this model at £ = 0 ,tt/T , 
yielding Wo = W n /T = 1 [24]. In a cylindrical geometry 
the two edges host linearly dispersing chiral modes in the 
quasi-energy gaps between the two bulk bands. 

We now introduce a specific form of a time-dependent 
disorder potential, U(£), which still allows for an exact 
solution. The full time-dependent Hamiltonian is given 
by Ho(t) = -Hclean (£) + V(t). During the fifth segment of 
the driving period, we let Vj (t) = JE Vj cj Cj , where Vj 
is a uniformly distributed in the range [— V, U], and Cj 
is the annihilation operator on site j. During segments 
1-4, V(t) = 0 . We choose V < 

By following the evolution of a state that is localized 
on a single bulk site j at time t = 0 , one can easily 
verify that this state is a Floquet eigenstate, whose quasi¬ 
energy is + Vj [here +1 (—1) refers to a site in the 
A ( B ) sublattice]. The Floquet spectrum consists of two 
bands, with quasi-energies in the range [±^ — V, + 
V]. One can similarly follow the evolution of a state 
that is initially localized on a site at the edge; in the 
geometry of Fig. 3 (viewed as a “strip” geometry with 
edges parallel to the horizontal axis), a state initialized 
on the A (B) sublattice at the top (bottom) edge moves 
by two lattice constants to the right (left) every driving 
period. Therefore, chiral edge states persist even in the 
presence of the disorder potential. 

As long as the gaps in the quasi-energy spectrum 
at £ = 0 ,tt/T remain open, the winding numbers at 
these gaps cannot change. Therefore, at least over a fi¬ 
nite range of the disordered potential strength, we have 
Wo = W n /T = 1 as in the clean limit [24]. Moreover, 
in the disordered system, all the bulk Floquet states 
are localized. Therefore, as argued in the previous sec¬ 
tion, the winding number is actually independent of the 
quasi-energy: W £ = 1 for all 5. We conclude that 
by the definition presented in Sec. II, the Hamiltonian 
Ho(t) = H c \ ean (t) + V(t) realizes the AFAI phase. 

Clearly, the above model utilizes a very specific form 
of the periodic driving and of the added disorder. Nev¬ 
ertheless, we argue that the AFAI is a robust phase that 
does not require fine-tuning. To demonstrate the robust¬ 
ness of the phase, we now consider a generic local per¬ 
turbation of Ho(t) that preserves the periodicity in time, 
H\(t) = Hoit) + AD(t), and show that the AFAI phase 
survives up to a finite value of A. 

The perturbation D(t) is assumed to be periodic in 
time and short-ranged in real space, such that the ma¬ 
trix elements of D(t) vanish beyond the rth neighbor on 
the square lattice. For V = 0 (no disorder) and A ^ 0, 


the bulk eigenstates of U(T) are generically dispersive 
and delocalized. However, we argue that for V > 0 and 
for a sufficiently small A, all the bulk Floquet states re¬ 
main localized. To see this, we derive a time-independent 
effective Hamiltonian H for the Floquet problem (on 
the torus) with V 7 ^ 0 ,A ^ 0 : 

e -iHfT = re -if 0 T dt[H 0 (t)+XD(t)]^ 

where T denotes time ordering. We further write the 
effective Hamiltonian as + D e ff, where 

corresponds to the unperturbed (A = 0) effective Hamil¬ 
tonian, defined such that its eigenstates he in the range 
[— ^). Here, we are considering a system with peri¬ 

odic boundary conditions; we will comment on the edge 
states later. 

The key point, which we show below, is that for suf¬ 
ficiently small V and A, the hopping matrix elements of 
the effective static Hamiltonian + D e ^ decay expo¬ 

nentially with distance. If, in addition, A <C V/£l, then 
all of the Floquet eigenstates remain localized [38]. 

To find the effective Hamiltonian for A 7 ^ 0 we need to 
solve for D e ff. The unperturbed effective Hamiltonian is 
of the form 


H 


eff _ 
(0) - 


E 


3 


(-1)^7r 
2 T 



c E’ 


( 12 ) 


where rjj =0(1) for j on the A(B) sublattice. We express 
D e ff as a power series in A, 

D eff = A£>$ + A 2 D™ + A 3 D&> + .... (13) 

To find we expand both sides of Eq. (11) in powers 
of A, and compare them order by order. The details of 
the calculation are given in Appendix B. The results can 
be summarized by considering the explicit representation 
of the operators as a tight-binding “Hamiltonian,” 

^?=i> (i4) 

h3 

Using the explicit form for given in Eq. (12), we find 
for the lowest-order term 


A (1) = 

ij e iE i:j T _ i 



(15) 


with V(t) = Uo(t,0)^D(t)Uo(t,0), where Uo(t,0) is 
the A = 0 evolution operator, and Eij = Vi — Vj + 
‘-jr 1 * is the zeroth-order quasi-energy differ¬ 
ence between the states localized at sites i and j. 

As long as EijT is smaller than 2tt for every pair of 
sites (which is the case for V < ^), the factor ~ 

in Eq. (15) is bounded. Similarly, the matrix elements 
are all non-singular (see Appendix B). Under these 
conditions, we expect the expansion in powers of A to 
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converge. In Appendix B, we argue that for sufficiently 
small A, the matrix elements of H decay exponentially 
with distance. Therefore, H has the form of a tight- 
binding model with random on-site potentials and weak, 
short-range hopping. In this context, we expect all states 
to remain localized up to a critical strength of A. 

Since all the bulk states remain localized as A is turned 
on, the chiral edge states that exist for A = 0 cannot 
disappear; the only way to remove them is by closing 
the mobility gap in the bulk, allowing the two counter- 
propagating states at the two opposite edges to backscat- 
ter into each other. Hence, we expect the edge chiral 
states, and the associated quantized pumping, to persist 
up to a critical value of A where the bulk mobility gap 
closes. 


V. NUMERICAL RESULTS 

Numerical simulations substantiate the conclusions of 
Sections II-IV. We will first briefly summarize our main 
findings, and then describe the simulations and results in 
more detail in the subsections below. For the simulations, 
a variant of the model discussed in Sec. IV, defined on a 
square lattice, is used: 

H(t) = H clean (t) + A D + J2 V i c h’ ( 16 ) 

3 

where H c \ edin (t) is the time-dependent, piecewise- 
constant Hamiltonian described in Sec. IV (pictured in 
Fig. 3). Using numerics, we are now able to study the 
more generic case in which the sublattice potential (de¬ 
note here by D), as well as the disorder potential are 
time-independent (in contrast to the model studied in 
Sec. IV). We define D = ^ and take Vj 

to be uniformly distributed in the interval [— V, V}. The 
parameters of the model are chosen to be A = 7r, and 
Sab = 0 . 

In the clean case (V = 0), the system exhibits an 
anomalous Floquet-Bloch band-structure: the Chern 
numbers of all the bulk bands are zero, but the winding 
number W £ = 1 for any value of 5 within each of the band 
gaps [24]. Such a band-structure is depicted in Fig. 2(b). 
When the disorder potential is turned on, however, the 
system enters the AFAI phase. Below, we show numeri¬ 
cally that the bulk states become localized, and coexist 
with edge states which occur in all quasi-energies. Fur¬ 
thermore, when the system is initialized with fermions 
filling all of the sites in the vicinity of one edge, while the 
rest remain empty, as in Sec. Ill, the disordered system 
exhibits quantized amount of charge pumped per period, 
when averaged over long times. Finally, we examine the 
behavior of the system as the strength of the disorder 
potential is increased. We find that when the disorder 
strength reaches a certain critical value, the system un¬ 
dergoes a topological phase transition where the winding 
number changes from 1 to 0. For stronger disorder, a 


“trivial” phase (where all bulk states are localized and 
there are no chiral edge states) is stabilized. 


A. Localization, edge modes, and quantized charge 
pumping in the AFAI 

The localization properties of the bulk Floquet eigen¬ 
states of (16) can be extracted from the statistics of 
the spacings between the quasi-energy levels. For lo¬ 
calized states, the distribution of the level-spacing is 
expected to have a Poissonian form. In contrast, ex¬ 
tended states exhibit level repulsion and obey Wigner- 
Dyson statistics [39]. To distinguish between these dis¬ 
tributions, it is convenient to use the ratio between 
the spacings of adjacent quasi-energies levels [40-42]. 
Choosing the quasi-energy zone to be between —tt/T 
and 7 t/T (i.e., choosing —i \oge l£T = sT for —tt/T < 
5 < 7 t/T), we label quasi-energies in ascending order. 
We then define the level-spacing ratio (LSR) as r = 
min {£ n , £ n +i} / max {£ n , £ n +i}, where 5 n = e n - £ n -i- 
This ratio, r < 1, converges to different values for ex¬ 
tended and localized states, depending on the symmetries 
of the system. For localized states, r\ oc « 0.39 [40], while 
for extended states, r ex t ~ 0.6 [42]. The latter value is 
obtained when one assumes that the quasi-energies are 
distributed according to the circular unitary ensemble 
(CUE) [42], and in the thermodynamic limit, coincides 
with the value obtained by the more familiar Gaussian 
unitary ensemble (GUE). 

Since the Floquet problem does not possess any generic 
symmetries such as time-reversal, particle-hole, or chi¬ 
ral symmetry, we expect its localization properties to be 
similar to those of the unitary class [43-45]. In analogy 
with the situation in static Hamiltonians in the unitary 
class [46], we expect that arbitrarily weak disorder is suf¬ 
ficient to localize the all Floquet states (on the torus). 
However, for weak disorder, the characteristic localiza¬ 
tion length £ can be extremely long, and easily exceeds 
the system sizes accessible in our numerical simulations. 
Therefore, the level spacing ratio is expected to show a 
gradual crossover from having the characteristic of delo¬ 
calized states, v ext ~ 0.6, when £ L, to the value that 
indicates localized behavior, r ext ~ 0.39, when £ <C L. 

This behavior is demonstrated in Figs. 4(a), panels (i)- 
(iii), where we plot the disorder averaged level spacing 
ratio r and the density of Floquet states, as a function 
of the quasi-energy for different disorder strengths. For 
weak disorder, VT = 0.5, panel (i) shows that the level 
spacing ratio is r ~ 0.6 in any spectral region where 
Floquet states exists. On the other hand, panel (iii) 
shows that already for VT = 4, the level-spacing ratio 
approaches r « 0.39 at all quasi energies, as expected 
from localized states. 

Note that, as the disorder strength increases, the level 
spacing ratio decreases uniformly throughout the spec¬ 
trum [Fig. 4(a), panels (i—iii); the same behavior is seen 
at weaker values of the disorder (not shown)]. There is no 






FIG. 4. Localization of Floquet states in the AFAI as a func¬ 
tion of disorder strength, computed for the model presented 
in Eq. (16). We use A = tt and an L x L system with peri¬ 
odic boundary conditions, (a) Quasienergy density of Floquet 
states per unit area (DOS) and level spacing ratio (LSR), for 
three values of disorder strength as indicated by the markers 
on the axis of panel (b). For all cases we take L = 70. (b) 
Finite size scaling of the localization transition. Level statis¬ 
tics in the delocalized regime are described by the Gaussian 
unitary ensemble (GUE), characterized by an average level 
spacing ratio r ex t ~ 0.60; in the localized regime, Poissonian 
level statistics give r\ oc ~ 0.39. These characteristic values 
are indicated by dashed lines. 


quasi-energy in which the LSR remains close to 0.6, cor¬ 
responding delocalized Floquet eigenstates. This is con¬ 
sistent with the expectation that the bulk Floquet states 
become localized even for weak disorder, and the local¬ 
ization length becomes shorter as the disorder strength 
increases. The behavior of the LSR as a function of sys¬ 
tem size, Fig. 4(b), also shows behaviour consistent with 
the above expectation. In contrast, if the bulk bands of 
the clean systems carried non-zero Chern numbers, delo¬ 
calized states would persist in the bands up to a critical 
strength of the disorder, at which point they would merge 
and annihilate. 

In the AFAI phase all the bulk states are localized, 
but the edge hosts chiral modes at any quasi-energy 
(cf. Sec. II). To test this, we simulate the time evo¬ 
lution of wavepackets initialized either in the bulk or 
near the edge of the system. We consider the system 
in a rectangular geometry. The initial state, IV’o)? is 
localized to a single site x 0 = (xo,yo)- To obtain in¬ 
formation on quasi-energy resolved propagation, we in¬ 
vestigate the disorder-averaged transmission probability, 

|Gw (x, xo, e) | 2 , which is a function both of quasi-energy 
e and the total time of evolution Tf = NT. Here, the 
bar denotes disorder averaging. The transmission am- 



Coordinate, x Coordinate, x 


(|Gjv(x, xo,e)| 2 ), (a.u.) 

FIG. 5. Wavepacket dynamics in the AFAI. Using the same 
model as in Fig. 4, we plot the amplitude of the transmission 
probability,^| Gw (x, xo, e) | 2 ), c.f. Eq. (17) obtained after a 
time-evolution of Tfi n = 300T and averaged over disorder re¬ 
alizations. We simulate a strip of size 20 x 100 with open 
boundary conditions, and plot (|Gw| 2 ) for several quasiener¬ 
gies e/Q = 0, T, |, (a) shows (|Gw| 2 ) when the initial 

wavepacket is chosen at the edge xo — (96,1). It indicates the 
presence of a robust edge mode at all the given quasi-energies, 
(b) shows the probability when the initial wavepacket is cho¬ 
sen in the bulk, x = (50,10). This indicates that the bulk 
Floquet states are localized. These simulations were carried 
out with a time step of dt — T/100. 

plitude in each disorder realization, Gw, is obtained by 
a partial Fourier transform of the real time amplitude, 
G (x, xo, t) = (x \U(t)\ V’o), and is given by 

1 N 

G n (x, x 0 , e) = — ^ G (x, x 0 , t = nT) e l£nT . (17) 

n =0 

The real time transmission amplitude G(t) is com¬ 
puted numerically by a split operator decomposition. 
Figs. 5(a),(b) show |Gw| 2 at different quasi-energies, for 
initial states on the edge and in the bulk, respectively. 
The simulations are done for a disorder strength VT = 4. 
At this disorder strength, the analysis of the level-spacing 
statistics shown in Fig. 4(a) indicates that all the bulk 
Floquet bands are localized with a localization length 
smaller than the system size. Fig. 5(a) shows the value 
of | Gw | 2 when the wavepacket is initialized at the edge 
of the system, x 0 = (1,1). The wave packet propagates 
chirally along the edge. The figure exemplifies that the 
edge modes are robust in the presence of disorder, and 
are present at all quasi-energies. Importantly, edge states 
are also observed at quasi-energies where the bulk den¬ 
sity of states is appreciable, indicating that the chiral 
edge states coexist with localized bulk states [the density 
of states in the bulk is shown in Fig. 4(a)]. 

In contrast, Fig. 5(b) shows |Gw| 2 for a wavepacket 
initialized in the middle of the system. The wavepacket 
remains localized at all quasi-energies, as expected if all 
bulk Floquet eigenstates are localized. This confirms that 
the model we study numerically indeed exhibits the ba- 
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FIG. 6. Quantized charge pumping in the AFAI. (a) Cumula¬ 
tive average of the pumped charge per cycle in the limit of long 
times, [c.f. Eq. (9)], as a function of disorder strength. 
For VT > 5, the localization length is sufficiently smaller than 
the system size, and Q ^ approaches unity. The inset shows 
the finite size scaling of for VT = 8. (b) Cumulative 

average of the pumped charge for N periods, ( Q) nt /N , as 
a function of N . The disorder strength used was VT = 8. 
The approach to the quantized value can be fit to a power 
law ( NT)~ V with v — 1.72, see the log-log plot shown in the 
inset. In both panels, we averaged the charge pumped across 
all the lines running parallel to the y direction of the cylinder 
(see Fig. 2) and over 100 disorder realizations. The system 
size used L x x L y = 50 x 50. 


sic properties of the AFAI phase: fully localized Floquet 
bulk states, coexisting with chiral edge states which exist 
at every quasi-energy. 

Next, we numerically demonstrate the quantized 
charge pumping property of the AFAI. Using the model 
described above, we numerically compute the value of 
Qoo given by Eq. (9) for a single value of the flux, 6 X = 0. 
When computing Q^, we averaged the charge pumped 
across all the lines running parallel to the y direction of 
the cylidner (see Fig. 2), as well as over 100 disorder re¬ 
alizations. In Fig. 6(a), we show the cumulative average 
of the pumped charge per cycle in the limit of long times, 
Q ^ [c.f. Eq. (9)] as a function of disorder strength. At 
weak disorder, when the localization length is smaller 
than the system size, Q ^ is clearly not quantized. How¬ 
ever, for disorder strength VT > 5, the value of Q^ 
quickly tends towards unity. This agrees with the re¬ 
sults presented in Fig. 4(a.iii), which indicate that at this 
disorder strength, the localization length is substantially 
smaller than L = 70. Finite size scaling demonstrating 
that Qoo indeed asymptotes to unity in the thermody¬ 
namic limit is presented in the inset of Fig. 6(a). 

The value of the cu mula tive average of the pumped 
charge over N periods, ( Q) nt /N [c.f. Eq. (7)] is plotted 
vs. N in Fig 6(b), demonstrating its approach to Q^ for 
large values of N (i.e., at long times). As in panel (a), we 
averaged over all the lines running parallel to the y direc¬ 
tion, and over 100 disor der realizations. We examine the 
asymptotic behavior of ( Q) NT and find a power law be¬ 
haviour of the form ( Q) NT = Qoo TcN~ v with v = 1.72, 
shown in the inset of panel (b). Note that for a sin¬ 
gle disorder realization and a single vertical cut, (Q) NT 
is expected to exhibit an oscillatory behaviour with an 


envelope which decays as 1 /TV, see Appendix Cl. This 
expectation is indeed confirmed by our numerical simu¬ 
lations, as we show in Appendix D. In contrast, Fig. 6(b) 
shows a power-law behaviour with a power larger than 
1 and no oscillations; this is clearly the result of averag¬ 
ing over the frequencies appearing in ( Q) NT /N for each 
disorder realization and vertical cut. The above results 
numerically confirm the discussion in Sec. Ill, and con¬ 
clude our numerical analysis of the AFAI phase. 


B. Strong Disorder Transition 

For sufficiently strong disorder, we expect the AFAI 
to give way to a topologically trivial localized phase in 
which the winding number vanishes. We now analyze the 
transition between the AFAI and this “trivial” phase. As 
explained above, the winding number W £ can only change 
if a delocalized state crosses through the quasi-energy 5 as 
disorder is added. In the AFAI phase all of the bulk states 
are already localized. How does the transition between 
the two phases occur? 

Clearly, at the transition, delocalized states must ap¬ 
pear in the quasi-energy spectrum. As disorder is in¬ 
creased, the delocalized states must sweep the full quasi¬ 
energy zone, changing the topological invariant W £ as 
they do so. The transition from the AFAI phase to the 
trivial phase can therefore occur through a range of dis¬ 
order strength V~ <V< V^ + , where V~ is the disorder 
strength at which the first delocalized state appears, and 
U c + is the disorder strength at which all Floquet states 
are again localized, and W £ = 0 for all e. Below we will 
support this scenario using numerical simulations, and 
furthermore provide evidence suggesting that the transi¬ 
tion is of the quantum Hall universality class. 

We study the same model used in Sec. V A and exam¬ 
ine the level-spacing ratio, r, as a function of disorder 
strength and quasi-energy. For this model, our simula¬ 
tions indicate V~ « V c + , within our resolution (limited 
by the system size). In Fig. 7(a), we plot r, averaged 
over disorder realizations and all quasi-energies. We see 
that at disorder strength V C T & 40 the level spacing ra¬ 
tio reaches r « 0.6, indicating delocalization. On either 
side of this point, r approaches 0.39 as the system size 
increases, which indicates localization. The peak in the 
value of r as a function of disorder gets sharper for larger 
system size, which is a signature of a critical point of 
this transition. In Fig 7(b), we show that at disorder 
strength V c , the LSR is independent of the quasi-energy 
with r « 0.6 (for disorder strengths close to V c , we also 
find that the LSR is independent of the quasi-energy, 
but with r < 0.6). This indicates that all of the Flo¬ 
quet states have a delocalized character at this disorder 
strength, which leads us to conclude that V c = V~ « U c + . 

At the critical point, V = V c , we expect the wavefunc- 
tions to have a fractal character [47]. This behavior is 
manifested in the distribution of the inverse participation 
ratio (IPR), P 2 = IVK r )| 4 - We study the distribution 
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FIG. 8. Floquet spectrum for Thouless’ quantized adiabatic 
charge pump, a) Quantized adiabatic pumping in a ID system 
is manifested in chiral Floquet bands that wind around the 
quasienergy Brillouin zone (right and left movers are shown 
in green and orange, respectively), b) Outside of the adia¬ 
batic limit, uj > 0, counterpropagating states hybridize, and 
all Floquet bands obtain trivial winding numbers; quantized 
pumping is destroyed. 


FIG. 7. Transition from the AFAI into a trivial phase at 
strong disorder, (a) Average level spacing ratio as a func¬ 
tion of disorder strength. On increasing disorder strength, a 
transition is observed between two localized phases with de¬ 
localized levels at VT « 40. Here, the level-spacing ratio has 
been averaged over all quasi energies, (b) Level spacing ra¬ 
tio as a function of quasienergy and its comparison with the 
DOS, indicating that the entire Floquet band is delocalized, 
(c) Effect of finite size of on the distribution of the participa¬ 
tion of ratio, P 2 , at a given disorder strength, VT = 40. The 
system sizes used for the simulations are L x x L y — 40 x 40, 
70 x 70, 100 x 100. The shape of the curve does not change, 
indicating a critical phase, (d) Scaling collapse of the three 
curves with D 2 = 1.3, where for a critical phase it is expected 
that (P 2 )~L- D2 . 


of the IPR, V (log P 2 ) 5 among all the Floquet eigenstates 
and averaged over disorder realizations. Fig. 7(c) shows 
the distribution for different system sizes. We note that 
the shapes of the distributions for different sizes are simi¬ 
lar, a signature of criticality. In two dimensions, the aver¬ 
age value of the IPR at a critical point is expected to scale 
like (P 2 ) ~ P -jC>2 , with D 2 < 2 [47]. Fig. 7(d) shows the 
scaling collapse of all the distributions. From the collapse 
we find the fractal dimension, D 2 = 1.3. The inset in this 
figure also shows a linear scaling log(P 2 ) ^ —D 2 \ogL. 
The critical exponent D 2 we find in our numerical simula¬ 
tions is close to the value found for the universality class 
of quantum Hall plateau transitions [47, 48], D 2 « 1.4, 
indicating that the transition from the AFAI to the trivial 
phase may belong to this universality class. This is nat¬ 
ural to expect, since, like the quantum Hall transition, in 
the transition out of the AFAI phase a delocalized state 
with a non-zero Chern number must “sweep” through ev¬ 
ery quasi-energy, to erase the chiral edge states. We ex¬ 
pect that the AFAI transition can be described in terms 
of quantum percolation in a disordered network model, 
similar to the Chalker-Coddington model for the plateau 
transitions [49]. We leave such investigations for future 
work. 


VI. DISCUSSION 

In this paper we have demonstrated the existence of 
a new non-equilibrium phase of matter: the anomalous 
Floquet-Anderson insulator. The phase emerges in the 
presence of time-periodic driving and disorder in a two- 
dimensional system, and features a unique combination 
of chiral edge states and a fully localized bulk. Such a 
situation cannot occur in non-driven systems, where the 
presence of chiral edge states necessarily implies the exis¬ 
tence of delocalized bulk states where the chiral branches 
of the spectrum can terminate. In a driven system, the 
periodicity of the quasienergy spectrum alleviates this 
constraint, allowing chiral states to “wrap around” the 
quasienergy zone and close on themselves. 

One of the key physical manifestations of the AFAI is a 
new type of non-adiabatic quantized pumping, which oc¬ 
curs when all states near one edge of the system are filled. 
It is interesting to compare this phenomenon with Thou¬ 
less’ quantized adiabatic pumping, described in Ref. [33]. 

The complementary relationship between pumping in 
the AFAI and the Thouless case is best revealed by 
first viewing the Thouless pump from the point of view 
of its Floquet spectrum. In Thouless’ one-dimensional 
pump, a periodic potential is deformed adiabatically such 
that in each time cycle a quantized amount of charge 
is pumped through the system. In the adiabatic limit, 
the quasi-energy spectrum of the pump exhibits one pair 
of counter propagating one-dimensional chiral Floquet- 
Bloch bands, which wrap around the quasi-energy Bril¬ 
louin [see Fig. 8 (a)]. The (nonzero) quasi-energy wind¬ 
ing number of each band gives the associated quantized 
pumped charge [3]. Importantly, for any finite cycle time 
the two counter-propagating states hybridize and destroy 
the perfect quantization of the charge pumped per cycle 
[Fig. 8 (b)]. 

In a strip geometry, the AFAI can be viewed as a quasi- 
one-dimensional system. As discussed in Sec. II, the sys¬ 
tem hosts chiral edge states that run in opposite direc- 
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tions on opposite edges. Furthermore, as shown by the 
spectral flow (see Fig. 2), these counterpropagating chi¬ 
ral modes cover the entire quasi-energy zone, analogous 
to the counter-propagating modes of the Thouless pump 
[Fig. 8(a)]. Crucially, however, the counterpropagating 
modes of the AFAI are spatially separated and therefore 
their coupling is exponentially suppressed : no adiabatic- 
ity restriction is needed to protect quantization. Thus 
quantized pumping at finite frequency can be achieved 
in the AFAI phase. 


How is the AFAI manifested in experiments? First, the 
localized bulk and chiral propagating edge states could 
be directly imaged, for example in cold atomic or optical 
setups. More naturally for a solid state electronic sys¬ 
tem, the pumping current could be monitored in a two 
terminal setup. Unlike the case of an adiabatic pump, 
where a quantized charge is pumped at zero source-drain 
bias, to observe quantized charge pumping in the AFAI 
the chiral propagating states of one edge of the system 
would need to be completely filled at one end of the sam¬ 
ple, and emptied at the opposite end. We speculate that 
this can be achieved using a large source-drain bias. A 
detailed analysis of such non-equilibrium transport in a 
two or multi-terminal setup, as well as an investigation 
of promising candidate systems, are important directions 
for future study. 


The implications of our results go beyond those specific 
to the class of systems studied in this paper. As a direct 
generalization of our results, one can consider construct¬ 
ing anomalous Floquet insulators in different dimensions 
and symmetry classes. Floquet-Bloch band structures 
which generalize those of Ref. [24] can serve as a starting 
point for constructing such anomalous periodically driven 
systems. Going beyond the single particle level, an im¬ 
portant challenge is to understand how the properties of 
the AFAI change in the presence of interactions. An ex¬ 
citing possibility is to obtain a topologically non-trivial 
steady state for an interacting, periodically driven system 
[14-16, 18, 19, 50]. The common wisdom dictates that 
a periodically system with dispersive modes is doomed 
to evolve into a highly random state which is essentially 
an infinite temperature state as far as any finite order 
correlation functions are concerned [42, 51-53]. Our re¬ 
sults on the single particle level demonstrate that it is 
possible to obtain a topological Floquet spectrum with 
no delocalized states away from the edges of the system. 
It is therefore possible that such periodically driven sys¬ 
tems can serve as a good starting point for constructing 
topologically non-trivial steady states for interacting, dis¬ 
order (many-body localized) periodically driven systems. 
What types of topological steady states can be obtained 
by this method, and what are their observable signatures, 
will be interesting subjects for future work. 
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Appendix A: Construction of the deformed 
evolution operator on the cylinder 


In this appendix we construct the deformed evolution 
operator on the cylinder, U £ (t) of Eq. (5), which is used 
to demonstrate the existence of delocalized edge states in 
the AFAI. The deformation is designed such that at t = T 
it interpolates smoothly between U(T) (corresponding to 
the cylinder ) in the vicinity of the edges of the cylinder, 
and 1 in the bulk of the cylinder. We first define the 
family of operators iF(s) = a (lh s)|r)(r|, where 


a(y,s) = < 


, (y-£ i) 
’(*2-H)’ 


A L y-£ 1 -y) 

’ (£2-£i) ’ 


y < h 
h<y<h 

h<y<L v -t 2 (Ai) 

Ly £2 ^ y ^ Ly 

y>L y -i 1 . 


Here, we choose £ <C G <C £2 An where £ is the bulk 
localization length. Analogously to Eq. (1), the family of 
deformed evolution operators is defined as 

tie (t, s) = U(t) exp [itF^HfFis)} . (A2) 

Here, is defined as in Eq. (1), i.e., H^ s = ^ log U(T), 
where U(T) is the evolution operator for a full period on 
the torus. The deformed evolution operator corresponds 
to U £ (t) = U £ (t, s = 1). 

Note that, strictly speaking, U £ {t) of Eq. (A2) is not 
precisely of the block-diagonal form of Eq. (5). It still has 
exponentially small but non-zero matrix elements con¬ 
necting the different blocks. However, a second deforma¬ 
tion can take U e {T) to the form appearing in Eq. (5). 
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Appendix B: Perturbative derivation of the effective 
Hamiltonian 


Here, we outline the details of the derivation used 
to demonstrate the perturbative stability of the AFAI 
phase. We first examine the effective static Hamiltonian 
defined in Eq. (11), expressed as a power series in A. We 
insert Eq. (13) into (11) and expand both sides in powers 
of A. For the left-hand side, we obtain 


-iT(H° f + D e ff) 


(Bl) 


= U$ 


1 — i 


^ dt[XD^(t) + X 2 D^(t) + ... 
fj' dtdt 1 (A Dg(t) + ..) (A D$(tj + ..) 


where U° = e~ itH °« and D^\t) = (U^D^U?. Note 
that the eigenvalues of are only defined modulo 

27ikj/T (where kj is an integer). The form of (t) 
depends on the choice of kj , while the evolution operator 
does not. To fix this ambiguity, we choose Ej to lie in 
the range [—7r/T, tt/T). 

The right-hand side of (11), expanded in powers of A, 
reads 


U(T) = U 0 (T, 0) 


- iX f dtV(t) 
Jo 

— A 2 f dt f dt'V{t)V{t') + 
Jo Jo 


(B2) 


where Uo(t,t f ) = Texp i f* #o(£)J and T>(t) = 

Uo(0,t)D(t)Uo(t,0). 

Equating (Bl) and (B2), and using e ~ lTHeff = 
Uo(T, 0), we find that 

f dtD^g (t) = [ dW(t), (B3) 

Jo Jo 


and likewise 

rT 


f dtD^ (t) = —i f dt f dt'V{t)V(t') 

Jo Jo Jo 

+ i [ T dt f dt'DW(t)D%(t'), (B4) 

Jo Jo 


and so forth. 

To find Dtff explicitly, we express them in the “tight- 
binding” form (14); inserting this form into (B3,B4), and 
using the fact that H contains only on-site potentials 
and no inter-site hopping, we arrive at 


[ dtD is\ t )=[ dt E' 

Jo Jo 


JEijt A (n) t 
l,j 


= E 


p iEijT _ -j 

— _-A- n) ctc- 


(B5) 


Equating this expression for n = 1 to the right hand side 
of Eq. (B3) gives Eq. (15). From (B5) we see that, as 
long as EijT < 2 tt for every pair of sites, is non¬ 
singular. As E^T —)> 27r, may diverge for all n, and 

the expansion in A fails. This reflects the fact that, in 
general, a Floquet operator whose eigenvalues are spread 
throughout the quasi-energy zone cannot be generated 
by a static, local Hamiltonian. 

From the form of the right-hand side of Eq. (B3), we 
can analyze the maximum range of the hopping matrix 
elements in . We denote the maximum range of the 
hopping matrix elements in D(t) by r, where r = 1 cor¬ 
responds to a nearest neighbor hop, r = 2 to second 
neighbors, and so on. Since the matrix elements of the 
unperturbed evolution operator J7o(t, 0) vanish beyond 
second neighbor sites on the square lattice, we find that 
contains matrix elements whose range is at most 
r + 4. Similarly, from Eq. (B4), A-^ vanishes beyond 

the 2 r + 6th neighbor, and more generally, A^ vanishes 
beyond range n(r + 2) + 2. Hence, the matrix elements 
of T e ff at range n(r + 2) + 2 contain the exponentially 
small factor A n . 


Appendix C: Quantized charge pumping and the 
winding number 


In this Appendix, we show that for the AFAI, a non¬ 
zero value for the winding number W £ implies quantized 
charge pumping on the edge of the system. As in Sec III, 
we take an initial state with all sites filled in a strip of 
width £ near one edge of the AFAI, [see Fig. 2(c)], and 
the rest to be empty. 

To calculate the time-dependence of the pumped 
charge, we begin by deriving an expression for the instan¬ 
taneous current flowing across a longitudinal cut through 
the cylinder, i.e., across a line parallel to the 7/-axis. The 
corresponding current operator is found by first allow¬ 
ing a flux 0 X to be threaded through the cylinder. Next 
we pick a gauge where the gauge (vector) potential is 
nonzero only on the links connecting sites with x = L x 
to sites with x m 0. The net current flowing across the 
cut between x = L x and x = 0 is then described by the 
operator I x (t) = dH(0 x ,t)/dQ x , where H(6 x ,t) is the 
Hamiltonian of the system in the presence of the flux 0 X . 
Here the tilde denotes the cylindrical geometry. 

Below, we first (Sec. C 1) show that when averaged 
over many periods, the charge pumped approaches 
a quantized value equal to the edge topological 
invariant n e d ge , expressed in terms of the spectral flow 
on one edge of the system [Eq. (10)]. We then show 
(Sec. C 2) that Qoo is in fact equal to the bulk topological 
invariant W £ , given by Eq. (2) 
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1. Quantized charge pumping from spectral flow 


2. Quantized charge pumping and the winding 
number 


We start from Eq.(7), which gives the charge pumped 
during the time integral 0 < t < r. The initial state, 
which is a single Slater determinant in terms of position 
eigenstates, is given by a superposition of Slater deter¬ 
minants in term of Floquet states 

i*> = £^IHi°>- ( C1 ) 

<s jes 

Inserting this into the expression for the pumped charge 
in the interval 0 < t < r, Eq. (7), we get 

< q)t = E A *s A s 'EE / dt^j(t )i \ipk(t )). 

s,s f jeskes' x 

(C2) 

The double sum in Eq. (C2) contains both diagonal and 
off diagonal terms for the contributions of single particle 
Floquet states. Denoting each of the contributions by 
Qjkj and using id t \^(t)) = H{t)\^(t))^ the different terms 
in (C2) can be written as 

Qjk = jT di<^.(i)| [do x {H\Mt))) - HdeJMt))} 

= J dt{(ipj(t)\de x id t \ipj(t)) 

+ id t ((ipj(t)\) d e Jipk(t))} 

= i [ dtd t {^j{t)\d ex \^ k {t)). (C3) 

Jo 


We now show that Qoo, the charged pumped over a 
period averaged over long times, is equal to the winding 
number W, appearing in Eq. (2). We show that in the 
limit of a large number of cycles, TV, the average charge 
pumped per cycle contains a quantized piece equal to the 
winding number, plus a small correction that decays with 
the averaging time at least as fast as 1/AT. 

Recall that the many-body initial state that we con¬ 
sider is a single Slater determinant with electrons pop¬ 
ulating all sites in the strip of width t near one edge 
of the cylinder. At time t , the expectation value of the 
current (I x {t)) is given by the sum of contributions from 
each of these single particle states, propagated forward 
in time with the evolution operator U = U(6 x ,t) for the 
system with the threaded flux 6 X . Defining a projector 
Vi that projects onto all sites within the strip of initially 
occupied sites, the current is given by 

(4 (t))=Tr{u\e x ,t) 9H ^ X 'J± U(O x ,t)V e y (C7) 

Using Eq. (C7), the total charge pumped in the first 
cycle, (Q) = J^dt (. I x (t )), is given by 

(Q) = J\ Tr { frt (0 X , t) dH ^ x ’ ~ U{6 X (C8) 

Rearranging and using the chain rule, Eq. (C8) becomes 


According to Floquet’s theorem, the Floquet states 
can be written as \^pj(t)) = e~ l£jt \<j)(t)), where | <j>j(t)) = 
\4>j(t + T)) is a periodic function. Substituting this into 
Eq. (C3) we obtain for the diagonal terms, 

< C4 > 

and for the off diagonal terms, 

Qjk = i (e <( ^- £fc)T - l) (^(O)I^J^(O)), j ± k. 

(C5) 

Clearly, when computing the average charge pumped over 
AT periods, (Q) NT /TV, the off-diagonal terms will give a 
contribution which decays as 1/AT, while the diagonal 
terms will give the contribution which does not decay 
with AT, 

lun^ (Q) nt /N = Qoo = T rij , (C6) 

where rij is the probability for the j th Floquet state to 
be occupied, rij = 1^)- Averaging over the flux 

values, we obtain Eq. (10), which is the result we set out 
to obtain in this subsection. 


pT 

(Q) = J dt Tr{p e C7+ (d 0x (HU) - Hde x tf) }. (C9) 


In the thermodynamic limit, the current (I X (0 X )) is 
expected to be insensitive to the value of the threaded 
flux. Thus replacing the pumped charge by its value 
averaged over all 6 X , and using id t U = HU, we obtain 

pOiTT pT 

(Q)= 2 nj d9x J dtTr{v e (Wde x d t U + d t U^de x U)}, 

(CIO) 

or equivalently, through integration by parts, 


r 2 tt 


^ ~ 2t r J 0 d ° x 


dt Tr {V e (d t U^d 6x U - de x U ] d t U)) 


(Cll) 

Inserting UW = 1 and using (d\W)U = —Wd\U in 
each of the terms in the above equation gives 


i r 2v 

(Q) = ^ l de x 


dtn{(uido x u) [(Wd t u),v e ]}, 

(C12) 


where we used Tr{Ve [A, B] } = Tr{A [ B , Vi]}. 

We now examine the cases for which the commutator 
in Eq. (C12) is non-zero. Denoting A = (WdtU), the 
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matrix element (r|[H, Vi\\r') is nonzero in the two cases 

<r|[I,P,]|r'> = -<r|i|r'), V e |r) = |r) ,V t |r') = 0 

(r\[A,V e ]\r') = (r\A\r'), V e \r)=0, VeW) = W). 

(C13) 

To set up a convenient means for enforcing the condi¬ 
tions above, we introduce an auxiliary gauge transforma¬ 
tion under which the single particle states on the sites 
|r) = | x,y) transform as 

|r) -)• |r) , y<i 

|r> -> e ie « |r), £ < y < L y . (C14) 

We denote the unitary operator that applies this gauge 
transformation as Gq . Because Eq. (C14) defines a pure 
gauge transformation, the pumped charge cannot depend 
on the value of 0 y . Therefore we are free to average (Q) 
over all such gauges. Using [Ge y ,Vt] = 0 and defining 

A(dy) = Gl y {&d t U)Ge y and B(0 V ) = Gl y (U^d 0x U)G 9v , 
we thus obtain 

p 2 ir r 2 n pT 

(C) = ^y dOyj d9 x J dtTi{B(0 y )[A(9 v ),V e ]}. 

(C15) 

Importantly, Eqs. (C13) and (C14) can be expressed as 
[A(0 v ),P e ]*'ido y A(0 y ), (C16) 

whereby the pumped charge becomes 


Therefore, in Eq. (C17) we can replace A and B with A 
and £>, giving 

(Q) = -^fdoj dtTr{(U^de x U)de y (U^d t U)}, 

(CIS) 

where for brevity we denoted dO x dO y = d@, and united 
the integrals under a single integral sign. 

Inserting UU^ = 1 between the two terms in the 
trace in the equation above, and again using (d\W)U = 
— Wd\U, we get 

(Q) = --^fdQf dtTr{(U'do K U)(uid ev dtU (C19) 

-(U^d 6v U)(U%U))}. 

Using 

u ] de v d t u = d t (u^de y u) - d t u^de v u 

= (U'd t U)(U'd 0y U) + d t (U'do y U), (C20) 

we get 

(Q) = ^fde£dtTr{(U%U)[{U% w U), (U'de y U)\ 
-{U ] de x U)dt{U^dg y U)). (C21) 

For the moment we focus on the second term in the 
above equation. Integrating by parts gives 


(Q) = -^2 C de v rde x [dtTv{B{O y )de y A(O y )). 

(C17) 

So far, we have expressed the average current using 
the evolution U on the cylinder. Here we aim to obtain 
a bulk-boundary correspondence, relating the pumped 
charge to the evolution operator on a torus , i.e., a ge¬ 
ometry without edges. We consider a completion of 
the cylinder to a torus with fluxes 0 X and 0 y threaded 
through the two holes of the torus. The torus Hamilto¬ 
nian, H = H(O x ,0 y ,t), is identical to G\ H(6 x ,t)Go y in 
the interior of the cylinder. The corresponding evolution 
operator is denoted by U = U(6 x ,6 y ,t). 

Importantly, V and A = Wd t U, as well as U and A = 
W dtU, are local operators for 0 < t < T. Therefore, up 
to corrections which are exponentially suppressed in the 
size of the system, ido y A = id$ A(6 y ), i.e., the derivative 
with respect to 0 y gives an identical result in the case of 
a torus and a cylinder. Moreover, since id$ A(Q y ) is also 
a local operator, the only matrix elements of (r'\B(O y )\r) 
contributing in Eq. (Cl7) are those for which r and r' 
are in the interior of the cylinder and close to the edge 
of the initially filled strip, i.e., y ~ £. For these matrix 
elements, B (defined on the cylinder) and B = U^de x U 
(defined on the torus) are identical (up to corrections 
which are exponentially small in the size of the system). 


dG dtTr\(U^de x U)dt(U^d ei 


= fdej dtTr{±d t ({U^de x U){U^de y U)) 

- ±(dtrfde m U)(U% y U) - ±(U'dtde m U)(U% y U) 

+ ±{U'de x U)(a t U% v U) + ±{U'do x U)(U'a t de v U)}. 

Using the cyclic property of the trace, WU = 1, and 
integration by parts with respect to 0 X and it is pos¬ 
sible to show that the third and fifth terms (containing 
the double derivatives) cancel. The second and fourth 
terms, using Ud\U^ = —d\UU^, can be shown to give 
an identical contribution to the first term in Eq. (C21), 
but with a factor of — Defining the functional W[U (£)] 
for a bulk evolution U(t) as 


W\U] = 


dtn{(u'dtU) [(u^d 0x u), (u'de v u )\}, 


the net charge pumped during one driving cycle, assum¬ 
ing initial filling of a strip of sites covering one edge, is 
given by 

(Q) = W [U] -f^J'dtd t Tr{(U% x U)(U'de v U)}. 

(C24) 
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It is important to note that W[U] is quantized (and 
equal to a winding number as discussed in Ref. [24]) only 
for the case where the evolution is periodic, satisfying 
U(T) = U( 0). For such “ideal evolutions,” the second 
term in Eq. (C24) clearly vanishes, and therefore the 
pumped charge is quantized and given by the winding 
number. 

For a “non-ideal evolution,” where U(T) ^ P(0), 
W[U] need not be an integer. However, if the initially- 
filled strip near the edge is wide enough such that all 
edge states are occupied with probabilities exponentially 
close to 1, then the spectral flow arguments presented 
in Sec. C 1 indicate that the average charge pumped per 
cycle will yield a quantized value, with a correction that 
vanishes at least as fast as 1/TV. As we now show, this 
behavior can be seen directly through further manipula¬ 
tions of Eq. (C24). 

Consider a “continued” evolution U(t), defined on a 
larger time period of 2T. We define U(t) such that 
it is equal to the original evolution operator U(t) for 
0 < t < T, and to e iH **^~ T ^U(T) = e^ eff ^“ 2T ) for 
T < t < 2T. As in Ref. [24], P eff = ^ log C7(T); in the 
discussion below the choice of the branch cut of the log 
is unimportant, as long as the system is in a localized 
phase. As constructed, U(t) is an “ideal” evolution in 
the larger period 2T, i.e., U(2T) = P(0) = 1. 

Starting with Eq. (C24), we add and subtract the 
quantity W[e lHett ( t ~ T '~\, i.e., Eq. (C23) with U(t) re¬ 
placed by e lHef{ ( t_T ). After a shift of the time variable by 
T, the added piece combines with W[U] to give W 2 [£7], 
where HA is defined as in Eq. (C23) with the time integra¬ 
tion taken from 0 to 2T. The subtracted piece remains as 
a correction. The charge pumped over one cycle is then 


(Q) = W 2 [U] - W[e 

- gh^deTr{(^^)(^^ ) }|^ (C25) 


where the last term arises from the full derivative (second 
term) in Eq. (C24). Crucially, because U is 2T-periodic, 
HA [(7] is a true winding number and is quantized. The 
issue remains to characterize the contributions of the sec¬ 
ond and third terms; below we show that they can be ne¬ 
glected in the limit of a large number of pumping cycles. 

Consider the average charge pumped over TV driv¬ 
ing cycles, = T J^ T dt (I x (t )). To analyze this 

quantity we repeat the manipulations leading up to 
Eq. (C24) above. In moving from the evolution oper¬ 
ator on the cylinder to that on a torus, see discussion 
above Eq. (C18), we furthermore use the fact that in the 
localized phase U and A = Ud t U remain local even at 
long times. Likewise, U and A = Ud t U (for the cylinder) 
are local in the y coordinate. In this way we find 


(Q) NT 

TV 


fw N [u] + 


m 

TV 


(C26) 


where HAv['] is defined in the same way as W\-\ in 
Eq. (C23), but with the time integration taken up to 


NT rather than T. The factor /(TV) on the right hand 
side of this equation arises from the term corresponding 
to the full derivative term in Eq. (C24); its magnitude is 
bounded, and therefore the ratio /(TV)/TV decays to zero 
as TV goes to infinity. 

To see how the last term in Eq. (C26) vanishes for 
large TV, consider U in terms of its spectral decomposi¬ 
tion, U(NT) = ^2 n e~ l£nNT P ni where P n is the projec¬ 
tor onto Floquet state n. Each derivative contributes two 
terms: U^dg^U =-i(^.£r n )7VTP„ + P n ^P„. We con- 
sider each of the four resulting terms from the product 
Tr[(Wdg x U)(Wde v U)\ separately. 

First, when both derivatives act on the quasienergies, 
we get (NT) 2 do x £ n de y e n . A nonzero value for these 
terms would imply the existence of a current that grows 
linearly in time, which is unphysical. Moreover, as shown 
in Appendix Cl, as a general rule the time-averaged 
current (or pumped charge) must limit to a constant 
plus a correction that decreases at least inversely with 
time. In the fully localized phase, the quasienergies { e n } 
are exponentially insensitive to changes in the fluxes 6 X 
and 0 y (see arguments below), and therefore these terms 
clearly give a vanishing contribution in the thermody¬ 
namic limit. Therefore these terms can (and must) be 
dropped within the level of all other approximations of 
exponential accuracy employed above. 

Next, when one of the derivatives acts on the quasi¬ 
energy and the other acts on a projector, we get terms 
like NT do x e n Tr [P n do y P n \- These terms strictly van¬ 

ish due to the general identity Tr[P^] = 0, for any pa¬ 
rameter A upon which the projector P depends. 

Finally, when both derivatives act on the projec¬ 
tors we get a nonvanishing contribution of the form 
En,m Tr [( P n9^P n )(P m 9^P m )]. Crucially, these terms 
do not depend on the length of the averaging interval, 
NT. Therefore the quantity /(TV) in Eq. (C26) is in fact 
constant in TV, and the ratio /(TV)/TV decays to zero in 
the long time (large TV) limit. Furthermore, in the local¬ 
ized phase, the contributions from projectors onto states 
localized far from the bonds where the gauge fields 0 X:V 
act are exponentially suppressed. Thus it is clear that 
for any fixed TV the quantity /(TV)/TV remains finite in 
the thermodynamic limit L x ,L y —> 00 . 

To evaluate HAv [U] in Eq. (C26) we break up the in¬ 
tegral over the range 0 < t < NT into TV segments of 
length T. Shifting the time variable within each segment 
to run between 0 and T, we obtain 


N—l 

W N [U] = J2 W [Un], U n {t) = U{t)U{nT). (C27) 

n =0 

As discussed for the “non-ideal” evolution U(t) above, 
the operators {U n } are not periodic in time and therefore 
VF[P n ] is not quantized. 

To isolate the quantized contribution to Eq. (C26) 
we add and subtract a “return map” contribution 
(nT)] for each term W[U n \. We further 
define the “continued” evolution U n (t) = U(t)U(nT ), 
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with U(t) as given above. Note that U n is periodic in time 
with period 2 T (though U n (0) = U n (2T) ^ 1), and there¬ 
fore W2[Un\ is separately quantized for each n. Moreover, 
by virtue of the fact that the winding number W 2 [U n ] is 
a topological invariant for periodic evolutions, its value 
cannot change under smooth deformations of U n . In par¬ 
ticular, we may deform U n —>> U via the continuous trans¬ 
formation U n (t ; s ) = U(t)U(( 1 — s)nT), by taking s from 
0 to 1. Hence we find that W 2 [U n ] = U^U], and there- 
foreE^WM =NW 2 [U]. 

Inserting the result above into Eq. (C26) and subtract¬ 
ing the appropriate return map contribution, we obtain 

Mp = W 2 [U] - 7 Wn [ e iH ^~ NT )] + (C28) 

where we have combined the contributions of the return 
maps for all n into one term Wn [ e zH effWNT)^ _ N 0 t e that 
by shifting time arguments we can make the replacement 
W N [e iH **«- nt)] = -W N [e -^*]. 

The quantity WN[e~ lHetft } is not necessarily quantized, 
since the unitary e~ lHeiit is not a periodic function of t 
over the range 0 < t < NT. However, if the eigenstates 
of i7 e ff are all localized, then we can show (see below) 
that WN[e~ lHeiit ] decays with N as 1/N (or faster). The 
underlying reason is that, in the localized case, the eigen¬ 
states of i7 e ff do not flow under insertion of the fluxes 0 X 
and 0 y into the torus. For now we simply assert this 
claim, and will prove it at the end of this section. Ac¬ 
cepting the claim to be true, we obtain 

Mp = w 2 [U] + /(AO /N, (C29) 

where f(N) is bounded by a constant as a function of N. 
Equation (C29) is the result we set out to prove in this 
Appendix. 

Finally, to close the loose ends, we show that 
Wn is bounded as a function of N. Using the 

spectral decomposition e ~ lHeiit = e~ lEnt P ni we have 

W N [e~ iHeiit ] 

_ r NT e -iAet £ r . 

= YJ dOdt — 2 w I 'r{P n [P mi de x P ma ,Pk 1 de v Pk,}\, 

(C30) 

with As = e mi + 5/ Cl — £ m2 — Sk 2 , and the sum taken over 
the integers n,mi 5 2 ,&i, 2 - To get to Eq. (C30), we used 

Tr|p ra [P m ,P fc ]| = 0 and Tr(P n [P m , P kl dg y P k2 ] } = 0. 
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When is fully localized, its eigenstates do not 
“wrap” around the cycles of the torus, and therefore are 
insusceptible to the flux insertion. Therefore, up to cor¬ 
rections which are suppressed as exp(—L/£), where £ is 
the localization length, we have that: (i) the eigenvalues 
£ n (0) are independent of the values of the fluxes; (ii) un¬ 
der changing the values for the fluxes, the projectors P n 
transform as if transforming under a local gauge trans¬ 
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Here, Q x and Q y are projectors on sites that define the 
gauge transformation felt by the localized eigenstates. 

Returning to Eq. (C30), we note that in order for 
a term in Eq. (C30) to grow with A, it must have 
As = 0. Excluding the possibility of a fine tuned de¬ 
generacy which occurs on a finite area in flux space, the 
condition A^ = 0 requires that either mi = m 2 , ki = k 2 
or mi = k 2 ,mi = k%. However, the contribution of 
the latter two cases to Wn [e ~ lHeiit ] can be shown to 
vanish by substituting do a Pm = ie zT [Q a , P m ] e _zT in 
Eq. (C30), and applying straightforward algebraic ma¬ 
nipulations. Therefore, we find that all the terms in 
Wn [e~ lHeift ] are bounded by a constant as a function 
of A, and thereby we obtain Eq. (C28). 


Appendix D: Charge pumping statistics 

As mentioned in the main text, quantization of the 
charge pumped per cycle is realized for every individual 
disorder realization in a large system. In this short ap¬ 
pendix we show numerical results for a single disorder 
realization in a finite system of size 50 x 50 lattice sites. 
In the main text we defined the pumped charge by in¬ 
tegrating the current across a single vertical “cut” across 
the system in a cylinder geometry (a cut along the y di¬ 
rection, c.f. Fig. 2). Here, in Fig. 9 we show that the 
detailed time dependence of (Q) NT i the cumulative av¬ 
erage charge pumped per cycle across each single cut, 
displays a unique pattern of decaying oscillations and 
limits to one at large times (panel a). Current conser¬ 
vation implies that the average currents across all cuts 
must be equal in the long time limit; the spread of values 
at large times is due to numerical discretization error. 
In Fig. 9(b) we show the current averaged over all ver¬ 
tical cuts in the same 50 x 50 system. Here the rapid 
convergence to the quantized value is clearly displayed. 
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FIG. 9. Quantized charge pumping in the AFAI in a single 
disorder realization for a system of size 50 x 50. All parameters 
are the same as in Fig. 6. (a) Cumulative average of pumped 
charge per cycle, Qoo, resolved for five different vertical cuts 
across the cylinder (cuts along the y direction, c.f. Fig. 2). 
The spread of asymptotic values arises from numerical error 
in the finite time differences used for evaluating the integral 
in Eq. (7). (b) The pumped charge, averaged over all longi¬ 
tudinal cuts in the same finite-sized system. For a thermody¬ 
namically large system, averaging over all cuts is equivalent 
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